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ABSTRACT 

A theoretical  model  is  developed  for  the  prediction  of 
ship  lateral  motions  in  oblique  seas.  The  asymptotic  behaviour 
of  this  model  in  waves  that  are  long  relative  to  ship  beam  is 
examined,  with  particul~r  emphasis  on  the  classical  problem  of 
rolling  in  beam  seas.  The  theoretical  prediction  of  roll  damp- 
ing is  discussed,  and  the  importance  of  including  dynamic  lift 
on  appendages  is  emphasized.  Fairly  extensive  comparisons  of 
predicted  and  measured  roll  response  are  made,  with  good  agree- 
ment at  all  headings  considered. 


SOMMAIRE 


On  a mis  au  point  un  modele  theorique  en  vue  de  pre- 
dire  les  mouvements  lateraux  de  navires  dans  des  mers  obliques. 
On  etudie  le  compor tement  asymptotique  de  ce  modele  dans  des 
vagues  qui  sont  longues  par  rapport  a la  largeur  du  navire,  en 
mettant  l'accent  par ticulierement  sur  le  probleme  classique  du 
roulis  en  mers  de  travers.  On  discute  de  la  prediction 
theorique  de  1 ' amor t issement  du  roulis  et  on  souligne  l'import- 
ance  d'inclure  une  balancine  dynamique  sur  les  appendices.  On 
a compare  un  assez  grand  nombre  de  resultats  en  ce  qui  a trait 
au  roulis  prevu  et  mesure  et  on  a remarque  un  bon  niveau  de 
concordance  a tous  les  caps  etudies. 
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INTRODUCTION 


Over  the  past  two  decades,  considerable  success  has 
been  achieved  in  the  theoretical  prediction  of  ship  heave  and 
pitch  motions.  Computer  programs  to  perform  such  predictions 
are  now  in  common  use,  and  a large  number  of  correlation 
studies  have  shown  that  the  accuracy  of  the  absolute  motion 
predictions  is  generally  very  satisfactory.  For  lateral  plane 
motions,  however,  of  which  by  far  the  most  important  is  roll, 
the  situation  is  much  different.  Although  several  programs 
exist  which  predict  roll,  sway  and  yaw  by  purely  theoretical 
means,  for  example  Reference  1,  correlation  studies^  have 
shown  that  errors  in  roll  prediction  are  generally  significant, 
particularly  for  high  speed  warship  hull  forms.  As  to  sway 
and  yaw,  correlation  studies  have  been  very  few  and  incon- 
clusive . 


The  fundamental  reason  for  the  discrepancies  reported 
in  correlation  studies  is  that  programs  such  as  Reference  1 
generally  make  inaccurate  estimates  of  roll  damping,  especial- 
ly at  higher  speeds.  Experience  with  these  programs  has  led 
to  the  widely  held  belief  that  the  roll  damping  prediction 
problem  is  so  complex  and  non-linear  as  to  be  theoretically 
intractable.  In  keeping  with  this  philosophy,  certain  com- 
puter programs^  require  that  the  roll  damping  coefficient  be 
input  by  the  user.  Model  test  programs  specifically  dedicated 
to  the  experimental  determination  of  roll  damping  are  fairly 
common . 


The  failure  of  Reference  1 to  make  accurate  predic- 
tions of  roll  damping  may  be  traced  to  an  inadequate  treatment 
of  hull  appendages,  and  in  particular  to  the  failure  to  in- 
clude the  effects  of  dynamic  lift  on  such  appendages  as  rud- 
ders, skegs  and  propeller  shaft  brackets.  The  theory  present- 
ed in  this  report  includes  these  effects  and  consequently  pro- 
duces good  estimates  of  both  roll  damping  and  roll  response. 
Demonstration  of  a conceptually  simple  yet  reliable  technique 
for  roll  damping  prediction  is  probably  the  single  most  im- 
portant contribution  of  the  present  report. 

The  theory  presented  herein  has  basically  four  facets: 

1)  Strip  theory  for  computing  hull  added  mass,  wave- 
making damping  and  exciting  forces. 

2)  Lifting  surface  contributions  to  damping  and  ex- 
ci ting  forces . 
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3)  Viscous  roll  damping,  principally  from  bilge 
keels  . 

4)  Hull  circulatory  effects. 

After  a detailed  presentation  of  theory,  asymptotic 
expressions  valid  for  long  waves  are  derived  and  discussed  in 
order  to  provide  a physical  "feel"  for  the  mathematical  model. 
Roll  damping  is  then  discussed  and  examples  are  presented 
showing  the  relative  significance  of  the  various  contributions 
to  roll  damping.  Finally,  theoretical  predictions  are  com- 
pared with  experimental  data,  principally  from  regular  wave 
tests  of  a fully-appended  self-propelled  frigate  model.  Ag- 
reement is  generally  good  and  it  is  concluded  the  theory 
yields  results  of  sufficient  accuracy  for  preliminary  design 
purposes . 


2.  EQUATIONS  OF  MOTION 

In  applying  strip  theory  to  a displacement  hull,  cer- 
tain assumptions  are  standard: 

i)  Ship  response  is  a linear  function  of  wave 
excitation . 

ii)  Ship  length  is  much  greater  than  either  beam  or 
draft. 

iii)  All  viscous  effects  other  than  roll  damping  are 
negligible . 

iv)  The  hull  does  not  develop  appreciable  planing  lift. 

The  linear  equations  of  motion  for  the  ship  without 
roll  stabilizers  are  given  below.  These  are  written  with  re- 
spect to  a stability  axis  system  fixed  in  the  ship;  this  axis 
system  is  slightly  different  from  the  translating  earth  axes 
employed  in  Reference  4.  Stability  axes  are  commonly  used  in 
dynamic  simulation  of  both  aircraft  and  marine  vehicles.  In 
the  present  case,  stability  axes  offer  the  advantage  of  being 
more  suitable  fo  control  studies  and  of  yielding  simpler  ex- 
pressions for  the  coefficients  of  the  equations  of  motion. 
Further,  hull  circulatory  effects  art  included  in  the  mathe- 
matical model,  and  expressions  given  in  the  literature  for 
these  effects  are  written  in  terms  of  stability  axes.  The 
axis  system  is  illustrated  in  Figure  1.  The  axes  are  fixed 
in  the  ship  with  origin  at  the  CG,  and  in  the  reference  con- 
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dition  of  steady  forward  speed  with  no  seaway  or  control  dis- 
turbances, the  x-axis  is  directed  horizontally  forward,  the  z- 
axis  vertically  upward,  and  the  y-axis  to  port. 


The  coupled  sway,  roll  and  yaw  equations  follow,  with 
notation  the  same  as  in  Reference  4. 


Sway:  (A^  + m)n2  + B^l^  + A,,^  + 

+ A26^6  + (B26  + mU)h6  = F 2 (1) 


Roll:  A24fi2  + B24n2  + (A44  + I4)n4  + B44h4  + c44n4 


+ A 4 6 ^ 6 + B46n6  = r4  (2) 


Yaw:  A62n2  + B62n2  + A^f,^  + B^r^ 


+ (A66  + I6)V,6  + B66\  = F6  (3) 


where  n is  sway,  r\ , roll  and  n,  yaw.  A.,  and  B.,  are  the 

2 4 b j k j k 

added  mass  and  damping  coefficients,  C are  the  restoring  co- 

J 

efficients,  and  F are  the  exciting  forces  and  moments.  These 
coefficients  are  ascribed  the  general  forms: 


. + 
I 


F C 

+ B . . + B . . 
iJ  ij 


where  superscript  H denotes  hull  coefficients  derived  from 
strip  theory,  superscript  F signifies  contributions  due  to  ap- 
pendages (foils)  such  as  rudders  and  fins,  and  superscript  C 
denotes  hull  circulatory  terms. 

For  B.,,  the  roll  damping  coefficient,  there  is  an  ad- 
4 A 

ditional  term,  B ^ , the  viscous  roll  damping  coefficient. 

This  takes  account  of  the  viscous  resistance  to  rolling  of 
bilge  keels,  skeg,  hull,  rudder  and  other  appendages. 

Detailed  expressions  for  the  terms  on  the  right  hand 
side  of  equations  (4)  to  (6)  are  given  in  Sections  3 to  5. 


D . HULL  COEFFICIENTS  DERIVED  FROM  STRIP  THEORY 

The  strip  theory  from  which  the  hull  coefficients  are 
derived  is  given  in  Reference  5. 


3.1  ADDED  MASS  AND  WAVE-MAKING  DAMPING 


The  expressions  given  below  for  hull  added  mass  and 
damping  are  considerably  simpler  than  the  corresponding  ex- 
pressions in  Reference  5.  This  results  from  using  stability 
axes  instead  of  translating  earth  axes.  The  transformations 
involved  are  discussed  in  Appendix  A. 
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The  above  integrations  are  over  the  length  of  the  ship 
In  practice,  the  hull  length  is  divided  into  approximately  20 
sections  and  the  two-dimensional  added  mass  (a.^)  and  wave- 


making damping  (bjk)  are  computed  for  each  section.  a2£ 


and 


b?2  result  from  swaying  motions,  a^and  apply  to  roll, 

while  a , and  b „ , are  due  to  cross-coupling  between  sway  and 
11  2 4 24 

roll. 

3.2  HYDROSTATIC  RESTORING  COEFFICIENT 


lateral 


The  only  hydrostatic 
mot  ions  is  C^^  : 


restoring  coefficient 


affecting 
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44 


mgGM 


(23) 


where  GM  is  metacentric  height. 


3.3  FORCING  FUNCTIONS 


The  strip  theory  of  Reference  5 yields  the  following 

expressions  for  the  sway  exciting  force  F? , the  roll  exciting 

H u 

moment  F, , and  the  yaw  exciting  moment  F" 

4 b . 


pa/ (f . + h . )d£  j = 2 , 4 

L J J 

(24) 

P a / [ ( f „ + h,H  + h9]d? 

T Z Z ICO  z 

(25) 

where  a is  the  amplitude  of  the  incident  wave  and  the  integr- 
ation is  over  the  hull  length.  f . and  h.  are  the  sectional 

J J 

incident  and  diffraction  forces,  respectively,  given  by 
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f.  is  also  commonly  referred  to  as  the  Fr oude-Kr i lof f 
force.  A physical  interpretation  of  f^  and  h^  is  that  f^  re- 
sults directly  from  the  action  of  the  incident  wave  system  on 
the  hull,  while  h^  represents  a correction  for  ship-wave 

interference  (diffraction).  The  integrations  are  performed 
over  the  submerged  hull  section  Cx>  n^  and  n^  are  the  y and  z 

components  of  the  unit  outward  normal  to  the  hull  at  (x,  y,  z), 
Figure  1 , and 


yn3  - zn2 


= z + h 


CG 


(28) 


(29) 


>2  and  <J>^  are  the  two-dimensional  section  potentials  for  sway 


and  roll  oscillations, 
number,  given  by 


to^  is  wave  frequency  and 


is  wave 


k = 
w 


(0 

w 

g 


(30) 


h is  the  height  of  the  centre  of  gravity  above  the  water- 

C G 

plane,  and  Bg  is  the  heading  angle  relative  to  the  sea, 
defined  in  Figure  2.  Note  that  for  head  seas  Bg  = 180  , for 
beam  seas  Bg  = 90°,  and  for  following  seas  Bg  = 0°.  to^  is 
related  to  to , the  frequency  of  encounter,  by 


to  = to 


k U cos  3 
w s 


(31) 


4.  VISCOUS  ROLL  DAMPING 


The  viscous  roll  damping  coefficient  may  be  expressed 

o f rt  1 1 nr.!o  • 


B44  = BBK  + BE  + BH  + BF 


where  B , B and  B denote  contributions  from  bilge  keels, 

B K.  E H 

eddy-making  resistance  of  the  hull  and  hull  friction,  respect- 
ively. B represents  the  viscous  effect  of  appendages  other 
r 

than  bilge  keels  (rudders,  fins,  etc.)  at  zero  speed.  Of  the 

four  components  of  bY.  , B is  normally  by  far  the  largest. 

4 4 dK 

Viscous  roll  damping  is,  of  course,  a non-linear 
function  of  roll  amplitude.  A linear  approximation  is  obtained 
by  equating  the  energy  dissipated  by  the  non-linear  viscous 
effect  during  one  roll  cycle  to  that  dissipated  by  a linear 
damping  term.  This  is  done  as  follows: 

Denote  by  B the  linear  damping  coefficient  and  let 
the  roll  angle  be  given  by 


H , = n , sinut 

4 4 


where  is  the  roll  amplitude.  Now  Br|^  represents  a torque 

about  the  CG  and  during  one  roll  cycle  the  following  amount 
of  work  must  be  done  against  this  torque: 


4/  ^BD/dr),  = 4a)2  b / 1I^'Jcos2ut  dt  = iruiBri 


If  the  energy  dissipated  by  the  viscous  effect  is 
E,  then  B is  given  by 


, 2 2 


BILGE  KEELS 


The  most  complete  study  to  date  of  the  efiect  of  bilge 
keels  on  ship  rolling  has  been  carried  out  by  Kato® , who  found 
that  the  following  factors  influenced  bilge  keel  e f f e c t i vne s s : 
bilge  keel  area,  breadth  and  aspect  ratio,  Reynolds  number. 
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ship  draught,  distance  from  the  bilge  keel  to  the  CG,  height 
of  the  CG  and  form  of  the  bilge.  After  analyzing  consider- 
able experimental  data,  he  devised  the  following  empirical 
method  for  calculating  the  energy  dissipated  through  bilge 
keel  act  ion . 

The  energy  dissipated  in  one  roll  cycle  is  given  by 


E 


4p£bkrn4 


rr>> 


C C C.C  (D  F 
o a k n w 


-ot 


(36) 


where  i is  the  bilge  keel  length,  b^  bilge  keel  breadth, 

r the  distance  from  the  center  of  the  bilge  keel  to  the  CG  and 

T is  the  period.  C , C , C.  , C , B and  F are  coefficients  de- 
r o a k ’ n 

pending  on  ship  form  and  Reynolds  number. 


A close  look  at  equation  (36)  shows  that  a viscous 
force  of  the  usual  form. 


Fv  • iDVRSBKCBK 


(37) 


acts  on  the  bilge  keel  with  moment  arm  r about  the  CG.  V 

R 

is  average  velocity  of  the  bilge  keel  due  to  roll,  S 

B K is 

bilge  keel  area  and  CDI,  is  a drag  coefficient.  Equations  to 

B K. 

compute  the  various  components  of  C (C  , C , etc.)  are  given 

B K.  O Si 

in  Appendix  B. 

From  equations  (35)  and  (36),  the  bilge  keel  damping 
coefficient  is  given  by 


BBK  ■ i>PtVS“VoCaCkC„®F~°‘  <38> 


In  evaluating  B the  same  approach  is  taken  as  in 
B K. 

computing  added-mass  and  wave-making  damping,  i.e.,  the 
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sectional  contributions  are  computed  individually  and  then 
s ummed . 


4 . 2 EDDY-MAKING 

The  results  of  this  section  are  due  to  Tanaka7,  who 
conducted  rolling  experiments  with  various  ship  sections  to 
assess  the  effect  of  section  shape  on  eddy-making  roll  damp- 
ing. 

The  r o 1 1- r e s i s t ing  force  due  to  eddy-making  is  ex- 
pressed as 


where  r is  the  distance  from  the  CG  to  the  point  on  the  hull 
where  the  eddies  are  generated,  S is  wetted  surface  area  of 
the  hull  section,  and  C is  a drag  coefficient  depending  on 
hull  form. 

Now  F exerts  a torque  about  the  CG  given  by 


The  energy  dissipated  by  this  torque  during  one  roll 


cycle  is 


, 1 a 4 , ~ , , 

E = 4 j Tdn^  = -jpr  SC 


whence,  from  (35),  the  eddy-making  damping  coefficient  is 


It  remains  to  evaluate  C.  Empirical  expressions  for 


7 = -rl  ( . 887  + .145C)(1.7T  + BCD)  + 2(KG  - T)  } 

D D 


where  C is  the  block  coefficient,  KG  height  of  the  CG  above 
B 

the  keel,  B ship  beam  and  T draft.  R^  is  a Reynolds  number 

based  on  average  rolling  velocity  and  average  distance  from 
the  CG,  7. 


4 . 4 APPENDAGES  OTHER  THAN  BILGE  KEELS 

At  zero  forward  speed,  viscous  drag  forces  opposing 
roll  act  on  appendages  such  as  rudders  and  fins.  By  regard- 
ing these  as  oscillating  flat  plates  and  applying  the  method 
described  by  equations  (33)  to  (35),  we  obtain  the  following 
expression  for  the  appendage  viscous  roll  damping  coefficient: 


BF  = l¥PWV(y2  + z2)3/2  SCn 


where  summation  is  over  all  foil  elements. 


is  the  normal- 


force  coefficient  for  a flat  plate  inclined  at  a large  angle 
to  the  flow;  Hoerner^  gives  a value  of  1.17  for  Cn  when  the 
angle  of  inclination  exceeds  40  . 


b.  LIFTING  SURFACE  CONTRIBUTIONS 

When  U > 0,  hull  appendages  such  as  the  rudder,  skeg 
and  propeller  shaft  brackets  act  as  lifting  surfaces  to  gener- 
ate both  damping  and  exciting  forces;  the  contribution  to  roll 
damping  is  especially  significant.  Of  lesser  importance  are 
the  appendage  added  mass  terms,  which  are  independent  of  speed. 

Since  derivation  of  the  lifting  surface  (foil)  co- 
efficients is  adequately  described  in  Reference  4,  only  the 
final  results  will  be  presented  herein.  Expressions  for  the 
foil  added  mass  and  damping  coefficients  are: 


r' 


' A 

» 

ft 

.1 

» 1 

•V 

*»■  I 


1, 


A = Ea  s in  2 T 

22  p 


B22  = |pu^scLac(k)sin2r 


A-2^  = -Ea  sinT(y  cosT  + z sin  f) 


B24  “ -^PUESC^^C (k ) s inF (y  cosP  + z sinT) 


26 


26 


Ea  s sin2r  = A^ 


62 


|pUESCLuC(k) (s  - c/4)sin2r 


A^  = Ea^  (y  cosT  + z sinF)  2 


~ |-pUESCLc(C(k)  (y  cosP  + z sinF)2 


= -Ea^s  sinf(y  cosT  4-  z sinf) 


(50) 


(51) 


(52) 


(53) 


(54) 


(55) 


(56) 


(57) 


(58) 


= -■|-pUESCL^C  ( k)  s inf  ( s - c/4)(y  cosT  + z sinF)  (59) 


62 


64 


^pUExSC,  C(k)sinzF 
2 Lu 


= — j pUE xS C ^ ( k ) s in T ( y cosF  + z sinF) 


(60) 


(61) 
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£a  s2sin2r 
P 


2pUISCLuC(k) (s 


- c/4)x  sin2r 


Summation  is  over  all  foil  elements,  including  rudder, 
skeg,  fins  and  propeller  shaft  brackets.  In  applying  these 
expressions  to  the  rudder  and  skeg,  the  dihedral  angle  F for 

fins  and  brackets  is  illustrated  in  Figure  3.  a is  the  added 

P 

mass  of  a foil  being  accelerated  perpendicular  to  its  surface. 
For  sufficiently  large  aspect  ratios,  such  as  the  A-brackets, 


a = TTpb(^)2 
P 2 


For  very  low  and  low  aspect  ratio  foils,  such  as  the  skeg  and 

rudder,  use  empirical  expressions  for  a . 

P 

A small  additional  contribution  to  B , . is  obtained  by 

4 4 

regarding  the  bilge  keels  as  very  low  aspect  ratio  lifting 
surfaces  : 


Bbr  = TipUb2  r 2 


where  b ^ is  bilge  keel  breadth  and  r the  distance  from  bilge 
keel  to  the  CG. 

The  exciting  forces  acting  on  the  appendages  are  ob- 
tained by  modifying  the  foil  exciting  forces  in  Reference  4 
to  apply  to  arbitrary  directions  to  the  sea. 


..  c F , P 

t ^sinf 


= - Z f „ ( y cosf  + z sinF) 
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(68) 


F r,F  . r 

F6  = slr*r 


where  f^  is  the  sway  exciting  force  acting  on  an  individual 
foil  element: 


-to  (sinFsinB  + i cosF)(— pUSCT  G(k)  + iwa  ) 
w s 2 La  p 


exp [-  k (h  + i(x  cosB  - y sinB  ))] 
w s s 


(69) 


The  methods  given  in  Reference  10  are  used  to  evalu- 
ate C(k),  the  circulation  delay  function,  and  G(k),  a modified 
form  of  Jones'  gust  function.  k is  the  reduced  frequency.  a 
is  foil  added  mass  and  is  included  in  equation  (69)  to  ac- 
count for  the  effect  of  wave  orbital  acceleration. 


6.  HULL  CIRCULATORY  EFFECTS 


Consider  now  the  hull  circulatory  terms.  These  are 
evaluated  following  the  methods  given  in  Reference  11.  Let 


B 


C 


(70) 


where  T is  draft.  Then,  from  Reference  11, 


B 


C 

22 


(71) 


B 


C 

26 


B„x 
C p 


(72) 


B 


C 

66 


B„(|c  L)2  + D/a  (d( 
C 2 p t 22 


(73) 
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where  x ts  the  x-coordinate  of  the  centre  of  area  of  the  hull 
P 

underwater  profile  and  is  the  prismatic  coefficient. 

The  action  on  the  hull  of  the  horizontal  component 
of  wave  orbital  velocity  generates  contributions  to  the  sway 
and  yaw  exciting  forces.  These  are  estimated  by  assuming  that 
the  hull  is  a wing  with  chord  L and  span  varying  with  section- 
al draft,  and  further  that  lift  is  uniformly  distributed  along 
the  ship's  length.  The  resulting  terms  are  given  below: 


f1'  = a T sinS  exp[-k  (ix  cos|?  + T / 2 ) ] (74) 

2 w x s w sx 

r B _ 

F2  = - ^ /f2dx  (75) 

P L 

g 

f£  = - ~ /f2x  dx  (76) 

P L 


7.  LATERAL  MOTIONS  IN  LONG  WAVES 


It  is  difficult  to  obtain  from  the  foregoing  mathe- 
matics a physical  appreciation  of  the  theoretical  model,  in 
particular,  the  expressions  for  the  hull  exciting  forces  are 
rather  abstract  in  appearance.  However,  when  wavelengths  are 
long  compared  with  the  ship's  beam,  fairly  simple  approxim- 
ations may  be  derived  for  the  hull  forcing  functions.  These 
asymptotic  expressions  provide  the  basis  for  further  theoret- 
ical analysis  of  the  lateral  motion  problem.  This  study  leads 
to  instructive  insight  into  the  relative  significance  of  the 
various  terms  in  the  equations  of  motion,  particularly  for  the 
classical  beam  sea  problem. 


7.1  ASYMPTOTIC  BEHAVIOUR  OF  THE  HULL  FORCING  FUNCTIONS 


Asymptotic  expressions  valid  for  waves  long  relative 
to  ship  beam  will  now  be  derived  for  the  hull  forcing  func- 
tions. Consider  first  the  sectional  wave  incident  force, 
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equation  (26) , and  expand  the  complex  exponential  in  the  in- 
tegral to  give 

I = /n.exp(k  z')(cos  p + i sin  p)d£  (77) 

C 3 
x 

where,  for  convenience. 


p = k y sinB 
r x s 


(78) 


From  symmetry  considerations, 
tegrates  to  zero  for  j = 2,  4 , since  n ^ 
board ) . Hence 


the 

(port) 


cos  p 
= -n 


term  in- 
(star- 


I = i/n.exp(k  z')sin  p d£ 
_ 1 w 

C J 

X 


(79) 


Now,  when  wavelength  is  long  relative  to  beam. 


sin  p = p,  exp (k  z')  = 1 
w 


(80) 


and  the  following  approximation  is  obtained 


I = ik_sinB„/n . yd£ 
'x 


w SC  3' 


(81) 


Now,  for  j = 2,  since  n^  = dz/d£ 


I = ik  sing  / y dz  = ik  sinB  A 
w s w s x 


(82) 


where  A is  sectional  area. 


For  j = 4,  from  equation  (28), 


I 


ik  sing  /( yn„ 
w sc  3 

x 


zn„)yd£  = -Ik  sing  M 
2 w s x 


(83 


where  is  the  sectional  contribution  to  the  r o 1 1- r e s t o r ing 

moment  (see  Appendix  E) . 

Combining  equations  (26),  (81)  and  (83)  yields  the 

following  long  wave  approximations  for  the  sectional  incident 
force  : 


f „ = -ik  g sing  A exp(-ik  x cosg  ) 
2 w s x w s 


(84 


f.  = ik  g sing  M exp(-ik  x cosg  ) 
4 w s x w s 


(85 


Consider  next  the  sectional  diffraction  force, 
equation  (27),  and  expand  the  complex  exponential  in  the  in- 
tegral. 


I 


f <b  . exp  (k  z ' 

C J 
x 


ftp  . exp  ( k z ' 

C 1 W 

x 


- n„smp  - 


) ( i n - - n>sing)(cosp+  i sin  p )d£ 
)[i(n^cosp  - n^sing^sinp) 
n sing^cos  p]  d £ 


(86 


From  symmetry  considerations,  the  term  in  brackets 
multiplying  i integrates  to  zero  for  j = 2,  4,  since 


j (port)  = -<p  . (starboard)  for  j = 2,  4 

n ^ (port)  = - n 2 (starboard) 


sinp  (port) 


-sin  p (starboard) 


Hence  , 


I = -/<(>. exp(k  z')(nsinp  + n.sing  cos  p )d£ 
C J w J z s 

X 


(87) 


In  long  waves,  the  approximations  (80)  apply,  as 

well  as 


cos  p = 1 


(88) 


.*.  I 


-f<p.  (n  k y sing 
C j 3 w 

X 


+ nosing  ) d£ 

z 9 


= -sing  /( i>  . n 0d £ 

sc  J 2 

X 


(89) 


since  k 

w 


is  small. 

By  definition^. 


I*Sn2i’-  * - 1“b2J> 

w 


(90) 


x 

Combining  (27)  and  (90)  yields  the  following  long 
wave  approximation  for  the  sectional  diffraction  force: 

-U> 


h . 
J 


— sin3  (b „ . + iwa  )exp(-ik  x cosB  ) 
3 s 2 j 2j  w s 


(91) 


7.2  LATERAL  MOTIONS  IN  BEAM  SEAS  OF  LONG  WAVELENGTH 


Specific  consideration  of  beam  seas  of  long  wave- 
length yields  considerable  insight  into  the  ship  lateral 
motion  problem.  This  is  important  from  both  a theoretical  and 
a practical  viewpoint  since 
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a.  classical  treatment  of  ship  lateral  motions 
(see,  for  example.  Chapter  IX  of  Reference  11) 
generally  concentrates  on  rolling  in  beam  seas; 

b.  for  ships  of  moderate  to  high  metacentric 
height,  rolling  is  greatest  in  beam  seas. 


Denote  by 
diffraction  forces. 


F_.  and  H.  the  total  wave  incident  and 
Ij  j 

respecively,  i.e.  the  itegrals  over  ship 


length  of  the  sectional  forces.  Now,  in  beam  seas,  equations 


(84),  (85)  and  (91)  become 


II 

CN 

4-1 

-ik  gA 

w X 

(92) 

f 4 = 

ik  gM 

w X 

(93) 

h . = 
1 

- >2J  + ‘"‘ij* 

(94) 

since  8 = 90° 

s 

and  a)  = a) . 
w 

Substitution  of  (92)  to  (94) 
evaluation  of  the  integrals  yields 

into  (24)  and  (25)  and 

FI2 

= -imgk  a 
w 

(95) 

H2  " 

-w(BK2  + i 0)  A ^ 2 ) a 

(96) 

FI4 

= ik  mgGMa 
w 

(97) 

Use  will  now  be  made  of  equations  (95)  to  (105)  to 
study  lateral  motions  in  beam  seas  of  long  wavelength,  be- 
ginning with  sway.  In  considering  the  sway  equation,  the  roll 
and  yaw  coupling  terms  will  be  ignored,  on  the  basis  of  the 
following  reasoning. 

The  roll-to-sway  coupling  terms  are  and 

note  that  in  common  with  all  wave-making  damping  terms, 

B £ ^ » the  dominant  component  of  B2^  . 0 as  u)  -*■  0.  Now  when 

waves  are  long  relative  to_ship  beam,  both  intuition  and  ex- 
perience show  that  n^,  = 0(£),  where 

w2  a 

C = -r-  (106) 

o 


The  argument  for  neglect  of  yaw-to-sway  coupling 
follows  similar  lines,  since  both  Intuition  and  experience 
suggest  that  yaw  Is  much  smaller  than  roll  in  beam  seas. 


Without  the  roll  and  yaw  terms  and  with  the  asymp- 


totic expressions  for  , equation  (1)  becomes 


-u)2(A22  + m)  + i0)Boo]no  = Fto  + H0  + F0  + F 


2 2 2 12 


a[-imgkw  - w(B22  + iwA22)  - U)(B22  + iwA22)  - ujB^  ] 


= ia  [-to 2 (A  „ + m)  + iwB00] 


22 


22 


(107) 


.'.  n2  = ia 


(108a) 


Thus,  sway  amplitude  equals  wave  amplitude  and  sway 
phase  leads  wave  phase  by  90°. 

Note  that  the  approximations  (95)  to  (105)  were  de- 
rived by  neglecting  the  exponential  decay  of  the  wave  forcing 
functions  with  depth  (equation  (80)).  If  this  effect  is  in- 
cluded, however,  by  simply  multiplying  the  right  hand  side  of 
equation  (107)  by  exp(-k  C T)  , where  T is  draft  and  Cw  midship 

W M M 

sectional  area  coefficient,  the  resulting  solution  for  r]  0 is 

0-  = ia  exp(-k  CWT)  (108b) 

2 w M 

which  is  valid  over  a wider  range  of  wavelengths  than  (108a). 

Figure  A demonstrates  that  equation  (108b)  gives  a 
very  good  approximation  of  sway  response,  especially  at  low 
wave  frequencies.  The  data  plotted  apply  to  a frigate  at 
cruise  speed  in  beam  seas.  The  solid  line  is  the  result  of 
solving  the  full  equations  of  motion,  while  the  dashed  line 
is  equation  (108b).  The  validity  of  neglecting  coupling  of 
roll  and  yaw  into  sway  is  clearly  demonstrated. 


Consider  next  roll  and  ignore  yaw  coupling.  From 
equation  (2)  with  the  asymptotic  value  (108)  for 


[_a)2(A44  + V + iwB44  + c44 

= FI4  + H4  + F4  ‘ ia(_a)2A24  + iwB24) 

= a[ikwmgGM  - u(b“4  + 1wa"4)  - w(BF4  + iaiA^)  + w(B24  + iajA24>  ] 

= ikwamgGM  = iC^?  (109) 


Thus,  sway  coupling  cancels  both  the  wave  diffract- 
ion moment  and  the  foil  exciting  moment,  and  the  result  is 
the  classical  uncoupled  rolling  equation.  Further,  as  id  4 0, 


n 


4 


U 


(110) 


Thus,  roll  amplitude  equals  wave  slope  and  roll 
phase  leads  wave  phase  by  90°. 

Consider  finally  the  yaw  equation.  From  equation 
(3)  with  n ^ given  by  (108), 


[-u)2(A..  + Ic)  + iwB66  ]n6  + [-w2a64  + iwB64]n4 


66 


H6  + F6  + F6  ' ia(-“2A62  + iwB62) 


= a|-a>(B26  + iuA|Jfi  + + 1u>a!J2  ) ] -u>BF6  - ioC 
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- i [-U)2  A„  , + UB»  + iw(B_,  + UA"  )]} 


Thus,  sway-to-yaw  coupling  cancels  all  the  yaw  ex- 
citing moments,  and  yaw  results  entirely  from  roll-to-yaw 
coupling : 


w2  A - iwB 

r 64 64 , 

n6  — 00 2 (A,  , + I.)  + iwB.,jn4 

bo  b b b 


(112) 


Since  the  denominator  in  equation  (112)  is  always 
much  greater  in  absolute  value  than  the  numerator  and 
H,  = 0(r, ),  H , is  very  small. 

4 D 

In  summary,  then,  the  foregoing  analysis  illustrates 
the  following  points  regarding  lateral  motions  in  beam  seas 
of  long  wavelength: 

1.  Coupling  of  roll  and  yaw  into  sway  is  of  little 
significance.  Sway  amplitude  approaches  wave 
amplitude  and  sway  phase  leads  wave  phase  by 
90°. 

2.  Sway  coupling  into  roll  cancels  the  wave  dif- 
fraction and  foil  exciting  moments,  and  the  roll 
equation  reduces  to  the  classical  roll  equation 
(109)  in  which  only  the  wave  incident  moment  ap- 
pears on  the  right  hand  side.  As  wave  frequency 
approaches  zero,  roll  approaches  wave  slope  and 
roll  phase  leads  wave  phase  by  90°. 


3.  Sway  coupling  into  yaw  cancels  the  yaw  exciting 
moment,  and  yaw  results  entirely  from  roll-to- 
yaw  coupling.  Yaw  is  consequently  virtually 
negligible  . 


4. 


The  derivation  of  equation  (109)  shows  very 
clearly  the  necessity  of  rigor  when  including 


the  foil  te  rms ; 

eluded  to  avoid 
frequencies . 


F F F 

A„, , B„.  and  F.  must  all  be  in- 

24  2 4 4 

erroneous  results  at  low  wave 


It  is  again  emphasized  that  the  above  considerations 
apply  only  to  beam  seas  of  long  wavelength.  In  particular, 
although  the  classical  roll  equation  provides  a good  approx- 
imation for  this  case,  it  is  wrong  to  assume  that  this  will 
be  true  at  other  headings  to  the  sea;  only  in  beam  seas  where 
wave  frequency  equals  frequency  of  encounter  does  the  can- 
cellation leading  to  equations  (109)  and  (112)  take  place. 


8. 


ROLL  DAMP  I N6 


Good  estimates  of  roll  damping  are  essential  for 
solving  the  roll  prediction  problem,  particularly  in  beam  seas 
where  for  long  wavelengths  equation  (109)  demonstrates  that 
rolling  dynamics  reduce  to  a simple  second  order  system.  It 
is  instructive  to  re-write  this  equation  in  the  conventional 

* 12 
form  : 


whe  r e 
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(114) 


is  roll  natural  frequency,  and 
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n = 


44 


oj  B . . 

o 44 


2 to  (A..  + I.) 
o 44  4 


2 C 


44 


(115) 


is  roll  decay  coefficient.  At  resonance  (oo  = ou  ) , equation 
(113)  gives 


n 


4 


2n 


(116) 


which  illustrates  the  importance  of  an  accurate  estimate  of 
n,  as  this  parameter  determines  peak  response  in  beam  seas. 

That  the  present  theory  does  give  reasonable  estim- 
ates of  roll  damping  is  demonstrated  in  Figures  5 and  6,  where 
theoretical  predictions  of  n are  compared  with  experimental 
values  for  two  different  warship  hull  forms  of  frigate  size. 
The  data  in  Figure  5 were  obtained  from  calm  water  experiments 
with  a 1:20. 8-scale  self-propelled  model  of  the  guided  missile 
1 3 

frigate,  FFG-7  , while  the  data  in  Figure  6 were  measured 

1 2 

during  full  scale  forced  rolling  trials  of  a Royal  Navy 
frigate.  Particulars  on  the  hull  forms  are  given  in  the  re- 
ferences. 
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In  Figure  5,  n is  plotted  against  initial  angle  of 
roll  for  speeds  of  0,  5,  10  and  15  knots.  The  measurements 
were  obtained  by  inclining  the  model  to  an  initial  roll  angle, 
then  releasing  it  and  measuring  the  resulting  roll  decay 
characteristics.  A large  number  of  measurements  were  taken, 
and  the  resulting  scatter  gives  a good  idea  of  the  experi- 
mental variability.  Theory  is  in  good  agreement  with  the 
measurements  at  all  four  speeds  shown.  Especially  noteworthy 
is  the  fact  that  the  rapid  increase  in  damping  with  increasing 
speed  is  well  predicted.  For  the  FFG-7,  by  far  the  greatest 
portion  of  this  damping  at  non-zero  speed  results  from  the 
large  spade  rudder  mounted  directly  behind  the  propeller.  In 
making  the  theoretical  estimates,  the  effects  of  propeller 
slipstream  were  accounted  for  using  the  data  presented  in 
Reference  14. 

The  ship  of  Figure  6 is  equipped  with  fin  stabiliz- 
ers, which  were  employed  in  the  trials  to  build  up  large  roll 
angles  which  were  then  allowed  to  decay  naturally  with  the 
fins  locked  in  the  mean  position.  The  resulting  measurements 
of  n are  plotted  against  speed,  together  with  the  theoretical 
estimates.  Again,  agreement  between  theory  and  experiment  is 
good  . 


The  importance  of  including  dynamic  lift  on  append- 
ages in  predicting  roll  damping  has  already  been  mentioned. 
Figures  7 and  8 illustrate  this  point,  giving  a theoretical 

V 

breakdown  of  B.,  between  viscous  (B  .)  and  non-viscous 
H F 

(B..  + B,.)  contributions.  Recall  from  Section  2 to  5 that 
4 4 4 4 

V F H 

®44’  ®44  ancJ  B44  are  » respectively,  the  viscous,  foil  dynamic 

lift,  and  hull  wave-making  components  of  roll  damping.  All 
damping  components  are  evaluated  at  the  roll  natural  frequency. 
Initial  roll  amplitude  is  5°  and  the  speed  range  is  0-25 
knots.  Figure  7 is  representative  of  a frigate  with  twin 
rudders,  one  pair  of  fins  and  one  set  of  bilge  keels,  while 
Figure  8 applies  to  a destroyer  with  a single  rudder,  no  fins 
and  a relatively  longer  pair  of  bilge  keels. 

The  dominance  of  the  non-viscous  contribution  at 
cruise  speeds  is  readily  apparent.  Further,  for  both  ships 
the  non-viscous  portion  results  almost  entirely  from  dynamic 
lift.  Figure  9 illustrates  this  point  for  the  destroyer, 

F H V 

giving  a plot  against  freqeuncy  of  B ^ , B^  and  B t ^ at  a speed 

of  20  knots.  At  the  roll  natural  frequency,  to  which  the  data 
in  Figure  8 apply,  wave-making  damping  is  relatively  low. 
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Figures  7 to  9 clearly  show  that  it  is  essential  to 
include  dynamic  lift  in  theoretical  estimates  of  roll  damping. 
Without  this  effect,  damping  will  be  significantly  under- 
predicted at  even  moderate  speeds. 


9. 


COMPARISON  OF  THEORY  WITH  EXPERIMENT 


Experimental  data  for  comparison  with  theory  have 

2 

been  obtained  from  three  sources:  regular  wave  tests  with  a 
fully-appended  self-propelled  frigate  model  at  the  seakeeping 
basin  of  the  David  Taylor  Naval  Ship  Research  and  Development 

Center  (DTNSRDC),  sea  trials^  of  the  oceanographic  research 

1 2 

ship  CFAV  QUEST,  and  sea  trials  of  the  same  ship  for  which 
roll  damping  data  are  given  in  Figure  6.  Only  roll  data  are 
presented,  since  measurements  of  sway  and  yaw  are  virtually 
non-exis  tent. 


Figures  10  to  14  present  extensive  experimental  data 
for  the  frigate  model  tested  at  DTNSRDC.  This  was  a compre- 
hensive experimental  program  in  which  metacentric  height, 
bilge  keel  configuration,  heading  and  forward  speed  were  sys- 
tematically varied.  Two  different  bilge  keels  were  used  (de- 


signated BK  and 


bk3) 


and  three  metacentric  heights  (GM 


1* 


GM, 


and  GM3) . 


Some  details  are  given  below. 


Length  between  perpendiculars 

Beam 

Draft 

Block  coefficient 
Prismatic  coefficient 

Metacentric  height  designation 
GM  1 

gm2 

GM3 

Bilge  keel 
designation 

BK  1 
BK. 


Area  in  % of 
wetted  area 

2 . 34 

1 .00 


18.182  ft 
2 .156  ft 
.716  ft 
.485 
.604 

% of  beam 
12  . 1 

9.0 

6 . 1 


Length  in 
% of  LBP 


29.6 

17.0 
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Figures  10  and  11  show  the  effect  on  roll  response 
of  varying  sea  direction  at  low  and  high  Froude  numbers,  re- 
spectively. Figures  12  and  13  show  the  effect  of  varying 
metacentric  height  and  bilge  keel  configuration  at  low  Froude 
number  in  bow  seas  and  high  Froude  number  in  beam  seas.  Fin- 
ally, Figure  14  presents  data  showing  the  influence  of  meta- 
centric height  on  roll  response  in  quartering  seas  at  low 
Froude  number . 

Agreement  between  theory  and  experiment  is  generally 
good  and  sometimes  near- pe r f e c t . The  following  points  are 
especially  worthy  of  mention: 

1.  The  considerable  variation  in  roll  response  with 
heading  to  the  sea  is  reasonably  well  predicted, 
both  at  low  and  high  speed. 

2.  Predictions  of  peak  response  are  generally  good, 
indicating  good  estimation  of  roll  damping. 

3.  Theory  correctly  predicts  a substantial  re- 
duction in  roll  response  as  speed  increases. 

This  decrease  is  mainly  due  to  the  twin  rudders, 
which  contribute  significantly  to  roll  damping 
at  high  speed. 

4.  The  dramatic  increase  in  roll  response  in  quar- 
tering seas  as  metacentric  height  is  reduced  is 
very  well  predicted.  This  is  of  considerable 
practical  importance,  since  the  G M ^ case  is  re- 
presentative of  the  metacentric  height/beam 
ratio  at  which  frigates  and  destroyers  commonly 
operate. 

Figures  15  and  16  present  sea  trials  data  for  CFAV 
QUEST,  for  which  the  leading  particulars  are^J: 

Displacement  during  trials  2255  tons 

Length  between  perpendiculars  234  ft 

Beam  42  ft 

Draf  t 16  ft 

Figure  15  compares  computed  and  measured  roll  re- 
sponse at  0 and  10  knots.  The  experimental  points  were  ob- 

tained from  roll  and  wave  spectra  measured  during  the  trials, 
simultaneously  in  the  case  of  zero  speed.  Agreement  between 
theory  and  measurement  is  good,  particularly  since  the  mea- 
surements come  from  full  scale  sea  trials  with  attendant  un- 
certainties in  wave  direction  and  spectral  measurement.  Again 
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the  decrease  in  peak  response  with  increasing  speed  is  cor- 
rectly predicted. 


1 


Figure  16  shows  that  when  the  theoretical  beam  sea 
roll  responses  given  in  Figure  15  are  applied  to  the  wave 
spectra  measured  during  the  trials,  a good  estimate  of  root 
mean  square  roll  angle  is  obtained.  Data  are  also  given  in 
Figure  16  for  the  bow  and  quartering  sea  directions,  and  again 
reasonable  agreement  is  demonstrated,  even  though  the  theoret- 
ical calculation  has  been  performed  for  a long-crested  sea. 

Finally,  Figure  17  presents  beam  sea  roll  response 
data  extracted  from  Reference  13.  Some  ship  particulars  are 
listed  below;  note  the  presence  of  a single  pair  of  fin 
stabilizers,  which  were  locked  in  the  mean  position  during 
the  trials. 


Displacement  during  trails 
Length  between  perpendiculars 
Beam 
Draft 

Stabilizer  fin  area 


2500  tons 
350  ft 
43  ft 
12  ft 
28  ft2 


Results  are  given  in  Figure  17  for  three  speeds  (12, 
16  and  20  knots).  As  for  the  QUEST  results  in  Figure  15,  the 
experimental  points  were  obtained  from  roll  and  wave  spectra 
measured  during  the  trials.  Three  sets  of  measurements  were 
obtained  at  16  and  20  knots,  and  two  at  12  knots.  Agreement 
between  theory  and  experiment  is  good,  particularly  at  the  two 
higher  speeds . 
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CONCLUSION 


A theoretical  model  has  been  developed  for  the  pre- 
diction of  ship  lateral  motions  in  oblique  seas.  The  asymp- 
totic behaviour  of  this  model  in  waves  that  are  long  relative 
to  ship  beam  has  been  examined,  with  particular  emphasis  on 
the  classical  problem  of  rolling  in  beam  seas.  It  has  been 
demonstrated  that  the  model  gives  reasonable  estimates  of  roll 
damping,  and  the  importance  of  including  dynamic  lift  in 
theoretical  predictions  of  roll  damping  has  been  emphasized. 
Fairly  extensive  comparisons  of  predicted  and  measured  roll 
response  show  generally  good  agreement  at  all  headings  to  the 
sea,  suggesting  that  the  theory  yields  predictions  of  suf- 
ficient accuracy  for  preliminary  design  purposes. 
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FIG.4,  FRIGATE  SWAY  RESPONSE  IN  BEAM  SEAS 
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FIG. 7,  FRIGATE  ROLL  DAMPING  COMPONENTS 
WITH  TWIN  RUDDERS, ONE  PAIR  OF  FINS 


FIG.8 , DESTROYER  ROLL  DAMPING  COMPONENTS 
WITH  SINGLE  RUDDER  , NO  FINS 
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FIG. 9,  DESTROYER  ROLL  DAMPING  COMPONENTS 

AT  20  KNOTS 
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FIG.  10,  FRIGATE  MODEL.HEADING  EFFECT,  Fn  = .I5,BK,  GM, 
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FIG.  1 1,  FRIGATE  MODEL, HEADING  EFFECT,  Fn  =.46, BK,  GM, 
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FIG.  12,  FRIGATE  MODEL,  GM  EFFECT,  Fn=.l5,BOW  SEAS 
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FIG.  14,  FRIGATE  MODEL, GM  EFFECT,  Fn  = .15 
QUARTERING  SEAS 
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FIG.  17,  FRIGATE  ROLL  RESPONSE  IN  BEAM  SEAS 
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FIG. 20,  FULL  SECTION 


FIG.2I. TRIANGULAR  SECTION 
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APPENDIX  A 


ADDED  MASS  AND  DAMPING  COEFFICIENTS  IN  STABILITY  AXES 

The  coefficients  given  in  Reference  4 apply  to 
translating  earth  axes.  The  translation  to  stability  axes  is 


n2  - n2  + un6  = n2  + ^ n6  - n2  - g2«6 


n4  * n4 


n6  - n6 


Consider  now  the  terms  on  the  LHS  of  the  sway,  roll, 
and  yaw  equations  in  Reference  4.  Ignore  the  end-effect  terms 
and  perform  the  above  axis  transformation. 


Sway  : 

[/a  dC](n2  + un6)  + [/b22dUOi2  - $rn6)  + A24n4  + B24n4 

L L 

+ tu22c*l  + £[b22dafi6  + '£b225d{  - U^22dt'\ 

= [/a  d51n2  + [/b22dUn2  + A24n4  + B24n4  + f/a22CdUn6 

L L L 


Roll: 


U 


[/a24df.)(n2  + Un6)  + [/b24d5](n2  - ^«6)  + *„n4  + »44fi4 

1-*  L 


u 


+ [/d24Cd?  + ^/b24d5)n6  + [/b24?d{  - u/a24dE]f,6 

L L L L 


[[a24dfJfi2  + llh2UA^2  + A 4 4 ^ 4 + B44%  + [T/a24^]H6 
L L L 


Yaw 


+ [/b24^]n6 

L 


u 


>fa225«  - ^b22dt>(d2  + ”V 

L L, 


+ r/b225d£  + u/a22dU(n2  - ^rfi6) 


+ A 6 4 ^ 4 + B64^4  + ^ {a22^2  d ^ + Z7!3  2 2' d U *6 

L<  L 


ir 


+ lfb22f  d«  + ^b22d«^6 

L»  L 


u 


• [,a22fdf-  - ^fb22df|d2  + ^b22ede  + ^‘22di^2 

L L L 


+ A 6 4 4 + B64*4  + lfa22^d^]fi6  + r/b22^d^6 

L Li 
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APPENDIX  B 


BILGE  KEEL  ROLL  DAMPING 

Equations  for  the  calculation  of  the  coefficients 

C , C m C,  , C , (§)  , F , and  a in  equation  (36)  are  given  in 

o a k n ^ 

this  Appendix. 

The  coefficient  F depends  upon  r,  , T,  and  T, 

the  angle  between  the  waterline,  CG  and  bilge  keel  root  (Fig- 
ure 18)  : 


F 


3 . 13rr). 

4 


rL7 


t/s\ 


The  index  a is  also  a function  of  r,  b,  and  T: 

k 


a = .6  - 2 . 03exp (-25C) 


whe  r e 


£ 


A_ 

rr  -75 


The  coefficient  (5)  depends  upon  the  length  of 
g rth  from  bilge  keel  root  to  waterline  s,  beam  B,  height  of 
CG  above  keel  KG,  draft  d and  rise  of  floor  F^  : 


® = cosy  + yf~[q  + P o - (p  q - p )f  (X)  ] 
k 
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where  y is  the  angle  made  by  the  plane  of  the  bilge  keel  with 
the  straight  line  passing  through  the  CG  and  the  bilge  keel 
root. 


r B / ^ £ \ , 

[-tan(^  - 2>  + Fr 


KG]sin(^-  + 


e = tan 


11  o ■ KG  - ! - K 


8 8 { KG  - d - . 54  [-^ 


(d  - Fr)tan(^  + y) ] } 


1 . 34sin 


f (X)  = 


ttX 

3.6 


1 + . 162sin  [j~g  (X  - . 9)  ] 


X 


2R) 


is  the  normal  pressure  coefficient  for  a rectang- 
ular plate  moving  with  a uniform  velocity  in  the  direction 
perpendicular  to  its  plane: 


C 

n 


1.  98exp(-llbk/«.  ) , bk/X,<.048 

1.17,  b.  Il>.  048 

k — 
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The  coefficient  depends  upon  ship  form  and  in 
particular  upon  R,  the  bilge  radius: 


C,  = 1 + 3 . 5e 
k 


■ 9k 


k = 


R ( i + X)2 


KG 


given  by 


Ca  is  a function  of  the  rolling  Reynolds  number  R^  , 


R 


N 


8bkrfU 

Tv 


where  V is  kinematic  viscosity. 


C 

a 


1. 


rn>103 


1.95  - . 251ogRN  + .2sin[-^(logRN  - 2.19)],  RN<103 


Finally,  the  coefficient  Cq  depends  upon  F and  also 
serves  to  scale  the  overall  equation: 


C = 14 . 1 + 37 . 3T3 
o 
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APPENDIX  C 


EDDY-MAKING  ROLL  DAMPING 


Empirical  equations  based  on  Tanaka's  results  are 
presented  below  for  evaluation  of  C,  the  eddy-making  drag  co- 
efficient . 


Consider  first  V or  U-shaped  sections.  Figure  19,  as 
are  encountered  in  the  foremost  part  of  a ship.  Tanaka  has 
obtained  the  following  empirical  equation. 


C = T .(zz)!  ,(a,y)exp(-uy)  (C-l) 

1 KG 


where  a is  the  angle  of  inclination  of  the  ship  side  at  the 
water  line  (Figure  19),  Rg  is  the  effective  radius  at  the  keel 

and  u is  a function  of  fi  . . T,  and  T„  are  tabulated  in  Tables 

4 1 2 

1 and  2 respectively.  Empirical  equations  for  Rg  and  u,  ob- 
tained by  fitting  least  squares  quadratics  to  Tanaka's  data, 
are  given  below. 


u = 14.1  - 46. 7n,  + 61. 7 ri  . 

4 4 


R = 
e 


|[4.12  - 3.69^  + . 823  (■— ■)  ] , ^<2.1 


(C-2) 


Consider  now  very  full,  almost  rectangular,  sections, 
as  shown  in  Figure  20.  Equations  (C-l)  and  (C-2)  are  used 
with  r the  distance  from  the  CG  to  the  bilge,  R the  bilge 
radius  and  T ^ = 1 . 


53 


Finally,  consider  triangular  sections,  Figure  21,  as 
are  found  farthest  aft  in  cruiser  stern  ships.  For  these 
sections  C is  a function  of  B/TCP  alone.  The  following  quad- 
ratic has  been  fitted  to  Tanaka's  data. 


C = . 438  - . 449  (^)  + . 236  (^;)2 

KG  KG 


B is  calculated  by  a strip  method,  with  each 
E 

section  either  placed  in  one  of  the  above  three  categories  or 
neglected.  An  example  of  a section  whose  eddy-making  effect 
is  negligible  is  a destroyer  mid-section  with  extremely 
rounded  bilges  and  no  skeg. 
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APPENDIX  D 


HULL  FRICTION  ROLL  DAMPING 

Consider  an  element  of  hull  surface  dS,  Figure  22, 

■+ 

The  skin  friction  drag  force  dF  acting  on  dS  as  a result  of 
the  rolling  velocity  f|^  is  given  by 


dF  = ^Pr(yn2  + z n_3  )n4l  n4f  CDFdS(n3T  _n2k) 


where  is  the  skin  friction  drag  coefficient,  j and  k are 

unit  vectors  along  the  y and  z axes,  and 


r = vy  + z‘ 


is  the  distance  from  dS  to  the  CG. 

dF  exerts  a torque  about  the  rolling  axis  given  by 


dT  = --pr(yn2  + z n3  )2  ^ | f)^  | CQ  pd  S 


The  energy  dissipated  by  dT  during  one  roll  cycle  is 


dE  = -pr(yn2  + z n^  ) 2 co2  fj  43CDpd  S 


Hence,  upon  integrating  dS  over  the  hull  by  a strip 
method,  equation  (43)  is  obtained. 
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APPENDIX  E 


ROLL-RESTORING  MOMENT 


Consider  an  element  of  hull  area  dA,  located  at  (x, 
y,  z).  When  the  ship  rolls  through  a small  angle  r\^,  the 
hydrostatic  force  acting  on  dA  increases  by 


df  = -pgn.ydA 


df  acts  normal  to  the  hull,  with  y and  z components: 


df  = -n.df 

y 2 


df  = -n  df 
z 3 

The  rolling  moment  exerted  by  df  is 


dK  = ydf z - zdf^  = (-yn^  + zn^Idf 

= -Pgn4y(-yn3  + zn2)dA 

= -pgr)4y(-yn3  + zn2)dxd£ 

Integration  over  the  hull  gives  the  total  roll- 
restoring moment 


K = -pgn,/M  dx 
4L  X 

where  is  the  sectional  contribution  to  the  r o 1 1- r es t or ing 

moment  : 

M -I  (~yn.  + zn  )yd£ 

X _ J £. 
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